Quantum melting of charge ice and non-Fermi-liquid behavior: 
An exact solution for the extended Falicov-Kimball model in the ice-rule limit 
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An exact solution is obtained for a model of itinerant electrons coupled to ice-rule variables on 
the tetrahedron Husimi cactus, an analogue of the Bethe lattice of corner-sharing tetrahedra. It 
reveals a quantum critical point with the emergence of non-Fermi-liquid behavior in melting of 
the "charge ice" insulator. The electronic structure is compared with the numerical results for the 
pyrochlore-lattice model to elucidate the physics of electron systems interacting with the tetrahedron 
ice rule. 
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The ice rule is a local constraint observed in a broad 
range of systems in condensed matter physics. It im- 
poses a configurational constraint on two-state variables 
denned at neighboring four lattice sites so that two out of 
four are in the opposite state to the other two. The most 
well-known material is water ice, in which the two states 
correspond to the configuration of hydrogens 0, [1J . An 
analogy was drawn by Anderson in the cation ordering 
of Fe 2+ and Fe 3+ in magnetite Fe3C>4 Q. More recently, 
a magnetic analogue was found in several pyrochlore ox- 
ides, the so-called spin ice, such as Ho2Ti 2 07 [4j and 
Dy 2 Ti 2 O r [|. 

The ice rule enforces local correlations; however, it is 
underconstraint and not enough to make the entire sys- 
tem ordered. The ground-state manifold retains macro- 
scopic degeneracy, resulting in residual entropy 0, 
Nevertheless, the ice-rule configuration is not completely 
disordered but cooperative in nature: There is a spa- 
tial power-law correlation in the two-state variables orig- 
inating from a hidden gauge structure @. Considerable 
progress on the understanding of such cooperative as- 
pects has been made in the last decade through the study 
of spin ice Q. 

In contrast to such "localized spin physics" , much less 
is known for the role of the ice rule in itinerant systems. It 
is intriguing to elucidate how the cooperative nature from 
the ice-rule constraint affects the electronic and trans- 
port properties. The issue will also be experimentally 
relevant to a wide range of pyrochlore-based compounds, 
such as mixed- valence compounds with a charge-ordering 
tendency [Iol - [l2j and itinerant (i-electron materials in- 
cluding Ising-like rare-earth moments [TH, Il4j . Only a 
few theoretical studies have been carried out so far Q . 

In this Letter, we address this issue in one of the sim- 
plest models which describe fermions interacting with ice- 
rule variables, an extended Falicov-Kimball model. We 
obtain an exact solution to this model on the Husimi 
cactus of tetrahedra, i.e., an analogue of the Bethe lat- 
tice composed of corner-sharing tetrahedra. We clar- 
ify the ground-state phase diagram including a "charge 
ice" insulator in which the fermions are localized in the 
ice-rule configuration. The solution reveals that a non- 
Fermi-liquid behavior appears at a quantum critical point 



(QCP) where the charge ice melts as the interaction de- 
creases. By comparison with the numerical results for the 
pyrochlore lattice, we show that our Husimi cactus model 
provides a cornerstone of itinerant ice-rule problems. 

We start with the extended Falicov-Kimball model on 
the pyrochlore lattice [Fig. QJa)], 

H = - < E( c |c 3 +H.c.) + C/^<(nf-i) + ^^nf^, 

where the first term describes the hopping of spinless 
fermions c, the second term represents the onsite repul- 
sion between spinless fermions and immobile particles 
/, and the last term is the repulsion between immobile 
particles. Here, = etc,, n{ = or 1 (c number), 




FIG. 1: (color online). A typical ice-rule configuration for (a) 
pyrochlore lattice and (b) tetrahedron Husimi cactus (THC). 
The sites with n{ = 1 (0) are shown by filled (open) circles. 
An example of loop and one-dimensional chain connecting 
the sites with n{ — 1 is shown by bold orange lines, (c) An 
apparently different ice-rule configuration obtained from (b) 
by interchanging the upper right and lower right branches, 
(d) A branch of THC considered in the calculations of g in 
Eq. ©. 
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and the sums are taken over the nearest-neighbor 

sites. Hereafter, we focus on the case in which the im- 
mobile particles satisfy the ice rule; see Fig.QJa). This is 
achieved by setting J^i n [ — N/2 (N is the total number 
of sites) and V — > oo. 

The partition function of the model is calculated by 
Z = Tr/Tr c exp(— f3H), where Tr/ (Tr c ) is the trace over 
the immobile-particle (spinless-fermion) degree of free- 
dom, and /3 — is the inverse temperature with the 
Boltzmann constant fee- For a given configuration of im- 
mobile particle the Hamiltonian ([1]) is reduced to 
a one-body problem given by 



n({n}) 



(2) 



where Ui is the onsite potential determined by the con- 
figuration {n{} as Ui — U(n{ — 1/2). Then the trace 
Tif in the partition function is replaced by the sum 

over {n{} which satisfies the ice-rule constraint: Z = 
S{nf } e ice Trc ^[-^({nf })]. This is, in principle, fea- 
sible to calculate because Tr c is performed by a diago- 
nalization of the one-body N x N Hamiltonian ([2]), but 
in practice, it is difficult for large system sizes because 
the sum within the ice manifold increases exponentially 
~ l.5 N / 2 i. 

A dramatic simplification is introduced by considering 
a modified structure of the pyrochlore lattice, that is, a 
Husimi cactus of tetrahedra [l5j]. It is an analogue of 
the Bethe lattice composed of corner-sharing tetrahedra, 
as shown in Fig. [lib), which we call here the tetrahe- 
dron Husimi cactus (THC). THC shares two important 
structural features with the pyrochlore — the tetrahc- 
dral units and their corner-sharing network. A difference 
is in the global connection of tetrahedra: The pyrochlore 
lattice has loops running across different tetrahedra [see 
Fig.QJa)], but THC does not have such global loops. De- 
spite of this difference, theTHC model captures several 
essential features of the pyrochlore, as we will see later. 

The simplification by considering THC is twofold. 
First, while the ice-rule configurations in THC are also 
macroscopically degenerate (= 6 x 3 N ^ 3 ), they are topo- 
logically equivalent because one can relabel the site num- 
bers by interchanging branches which spread from the 
same tetrahedron [see Figs.QJb) and (c)]. Consequently, 
all the possible ice-rule configurations give an identical 
Boltzmann weight, and therefore, the sum over {n{} can 
be suppressed in the calculations of any observable. 

On top of that, the second crucial point is that for 
any ice-rule configuration on THC we can obtain one- 
body Green's functions exactly by using recursion equa- 
tions similar to those often used in the Bethe lattice 
problems [l6l - [l8j . To see this explicitly, let us consider 
the T = local retarded Green's function at the site i, 
Gi(e) = t{i\[e — "H({n{}) + i<5]~ 1 |i) (t is included to make 
G dimensionless) . Similarly, we define the Green's func- 
tion gi (e) for a branch given by terminating a half of the 
tetrahedral network at the site i, as shown in Fig. [TJd). 



go is formally written by the expansion in terms of the 
hopping: 



(o) , (0)/ , , \ 
.9o = So + % (.9i + .92 + .93)50 

- 5o° } [51(32 + 33) + 32(33 + 3i) + .93(31 + 32)]3o 



(3) 



where <7g°^ = t(e — Uq)^ 1 is the atomic Green's function. 
Here, the second term corresponds to the processes where 
a fermion hops from the site i = to one of the other 
three sites in the same tetrahedron, j — 1, 2, or 3, then 
propagates within the branch belonging to the site j, 
and returns to the original site i — [see Fig. Did)]. 
The next term describes higher-order contributions, e.g., 
a hopping process — > 1 — > 2 — >• with propagations 
within each branch connected with the sites 1 and 2. 
The expansion $Ji§ is simplified and reduced to a set of 
recursion equations owing to the self-similarity of THC 
under the ice rule: gi depends not explicitly on i but 
only on the value of Ui] i.e., gi — g± corresponding to 
Ui = ±U/2. The recursive equations are given by 



3±(1-3t) + 2 3t( 1 "3±) 



l+.g T (l-2. 9± ) 
The full Green's functions G± are similarly obtained as 



G- ± l = 



9± 



I, U, 
^2) 



(5) 



Equations (Q} and ([5]) give the exact solution to the local 
Green's functions of the extended Falicov-Kimball model 



([T]) in the ice-rule limit: ^ 



N/2 and V -> 00. 



By extending the above calculations, it is also possible 
to obtain the nonlocal as well as the finite-T Green's 
functions. The whole procedure can be straightforwardly 
applied to a broader range of models with general ice-rule 
type constraints on general cacti of a complete graph. 
Such extensions will be discussed elsewhere. 

Now let us discuss the ground-state properties derived 
from the exact solution. Figure [2] shows the site-resolved 
density of states (DOS) given by p± = —lmG±/Trt and 
their summation p = + In Fig.[2ja), we show DOS 
at U — 0, which is given by the analytic form of p(e) = 
iRe{[ v /12< 2 -(e + 2t) 2 ]/[16t 2 - (e + 2t) 2 }} + \5{e - 2t). 
This is composed of the delta-functional peak (flat bands) 
at e = 2t and the broad spectrum for \e + 2t\ < 2 v / 3i- In 
the opposite limit of U 3> t also, a simple analytic form 
is available: By approximating g± ~ ±0 for e ~ TU/2 in 
Eq. fllj), we obtain p±(e) = iRe{4£ 2 - [s=p ([//2)] 2 } -1 / 2 . 
This is identical to DOS for one-dimensional (ID) tight- 
binding chains centered at e = ±U/2. The coincidence 
comes from the fact that THC under the ice-rule con- 
straint is broken up into ID chains of the same potential 
+U/2 or —17/2, as schematically shown in Fig. [Ub). A 
typical result is shown in Fig. [2te) at U = 10CW. 

In both the limiting cases of U = and U 3> t, the 
system is insulating at half filling of the spinless fermion, 



3 



n c = (nf) = 1/2. However, the origins of the two in- 
sulating states are quite different. The insulating phase 
at U = is a simple band insulator, in which the broad 
spectrum is fully occupied and the flat band is empty. On 
the other hand, the large-C/ phase is an incompressible 
state with a large gap ~ U originating from the repul- 
sive interaction between fermions and immobile particles. 
We call this correlation-driven insulating state "charge 
ice" , in analogy with the spin ice [3, @ , since the mobile 
fermions are localized in the ice-rule configuration which 
is composed of the sites with n( = 0. 

The question is how the system changes from the band 
insulator to the "charge ice" as U increases. In Figs.|S£b)- 
(d), we show the change of DOS obtained in the inter- 
mediate range of U. When switching on U, while the 
broad spectrum is almost unchanged, the fiat bands are 
perturbed to be broadened around e = 24, leading to a 
reduction of the gap [Pig. [2fb)]. The gap decreases as U 
increases and vanishes at U ~ 24 [Fig. d[c)]. Further in- 
crease of U opens a gap again [Fig.^d)], and increases it 
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continuously. Critical behavior around U = 24 is shown 
in Fig. Ufa). The plot suggests that the gap closes only 
atU = U C = 24. 

The critical behaviors of DOS can be understood an- 
alytically from Eqs. (j4| and (JSJ. We can prove that the 
divergence occurs at e — Srj = t+(U c — U)/2: At e = Etj, 
we obtain^ 1 = -l-(U-U c )/t, gZ 1 = 1, G+ = -t/U, 
and G- = oo, resulting in p + (su) = (except for U = 0) 
while p-(eu) — > oo. With regard to the gap, by con- 
sidering a small deviation from U = U c and evaluating 
g± at s ~ sjj, we find that the energy gap opens as 
A(U) ~ ^wlU - U c \ 3 for both U > U c and U < U c . 
Therefore the energy gap closes only at U = U c . This 
is identified as QCP between the band insulator and the 
charge ice (see the phase diagram in Fig. 01. 

The critical behavior of the gap A({7) oc \U — U c \ 3 is 
peculiar in contrast to the usual linear behavior A(C7) oc 
\U - U c \ in the Mott transition [20]. Actually, QCP is 
peculiar also in the sense that the self-energy exhibits 
an anomalous power-law behavior. From Eqs. Q and 
(O, we can derive that the self-energy E±, defined by 

,(0)-l_(J-l 



■9± 



shows the following critical behavior: 
3sgn(e - (i c ), 



ReE±(e) = 2-C±\e-(i, 
lcaS±(e) — — C±V3\e — p c 



(6) 
(7) 



where C_ = (4/t)s, C+ = C_/2, and fi c = t is the crit- 
ical chemical potential [Figs. O^b) and (c)]. The anoma- 
lous power law oc \e — /x c | 3 indicates that the system 
shows a non-Fermi-liquid behavior at QCP. DOS also 
shows a singular energy dependence, p±(e) oc \e — /x c | ±3 
at U — U c , resulting in the anomalous T dependence 
of thermodynamic quantities at QCP. For example, the 
specific heat is predicted to behave as oc T3 at low T. 

Collecting the results with varying U, we summarize 
the exact ground-state phase diagram in Fig. [4] There 
are two metallic regions for < n c < 1/2 and 1/2 < n c < 
1, which are separated by the two insulating regions at 
half filling n c = 1/2, i.e., the band insulator for U < U c 
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FIG. 2: (color online). DOS for (a)-(e) THC and (f)-(j) py- 
rochlore lattice at U/t = 0, 1, 2, 3, and 100. The results for 
the pyrochlore model are calculated for a 3 3 superlattice of 
4 x 8 3 lattice sites [ljj. p+(s) [p-(e)] is shown by the bold 
red (dotted blue) curves. In (a)-(d) and (f)-(i), the total DOS 
p(e) is plotted by thin black curves. 



FIG. 3: (color online). (a) DOS around e — p c — t 
at U/t = 1.6 (dash-double-dotted), 1.8 (dash-dotted), 2.0 
(solid), 2.2 (dashed), and 2.4 (dotted), (b) and (c) Non-Fermi- 
liquid behaviors of ImE± [Eq. J7J]. The green dashed lines 
show the asymptotic power law oc |e — p c \ 3 . 
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FIG. 4: (color online). Ground-state phase diagram of the 
extended Falicov-Kimball model on THC in the ice-rule limit. 



and the charge ice for U > U c . All these four phases 
meet at QCP at (U c ,(j, c ) = (2t,t). 

As indicated by the bold curve crossing QCP in Fig.[U 
the divergence of p„ at e = sjj is transferred from the 
upper-band bottom to the lower-band top [see Fig.^a)]. 
The divergence contains a bunch of extended states on 
the sites with n{ = 0. The extended states are fully 
occupied in the charge ice state, while they are empty 
in the band insulator. This fact leads us to define an 
"order parameter" to distinguish two insulating states, 
= (c\,cy) for the corresponding extended-state op- 
erator = Cj(— \yj\fL, where i represents the se- 
quential site number on a chain composed of n( = sites 
(L is the length). changes from for the band insula- 
tor to 1 for the charge ice. The transfer of the extended 
states bears some analogy to the "levitation scenario" 
proposed for the quantum Hall systems [2l|. This anal- 
ogy suggests a discrete change of the "transport nature" 
at QCP. Further analysis will be discussed elsewhere. 

Finally, we return to consider the original pyrochlore 
model. As indicated in Fig. [5J we observe many similar 



behaviors in DOS between THC and pyrochlore mod- 
els [l9|, [22j : (i) DOS consists of the fiat bands and dis- 
persive bands at U = 0, (ii) U drives a quantum phase 
transition to the charge ice insulator at half filling, (iii) 
the divergence of DOS transfers through QCP, and (iv) 
DOS shows a one-dimensional-like form for U ^ t. These 
are direct consequences from the key features shared be- 
tween THC and pyrochlore, i.e., the corner-sharing net- 
work of tetrahedra and the resulting macroscopic ice-rule 
degeneracy. Interestingly enough, (iii) suggests a possi- 
bility that the transition in the pyrochlore case is also 
described by an order parameter analogous to 0*, ac- 
companied by similar anomalous critical behavior. We 
note that, in the weak U region, THC is a band insu- 
lator, whereas the pyrochlore model is metallic at half 
filling [Figs. |2{f)-(h)]; however, this difference is rather 
irrelevant since our focus is on the strongly correlated 
physics related to the charge ice insulator. The bene- 
fit from obtaining the exact solution exceeds the minor 
dissimilarity. Thus, our THC solution captures the com- 
mon essential physics of the itinerant ice-rule systems. 
Further comparisons, including the effect of global loops 
neglected in THC, will be reported separately [22| . 

In summary, we have exactly solved the extended 
Falicov-Kimball model on the tetrahedron Husimi cac- 
tus in the ice-rule limit. The solution reveals a quantum 
critical point and associated non- Fermi-liquid behavior in 
quantum melting of "charge ice" insulator. Furthermore, 
the results capture many essential features of more real- 
istic lattices with corner-sharing tetrahedra, such as the 
pyrochlore lattice. Our exact solution provides a canon- 
ical reference to the itinerant ice-rule physics, and will 
open new avenues of research with wide applicability to 
correlation-induced phenomena under strong frustration. 
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